m  7^4  5578 


SAM-TR-72-10 


C9MPUTEit  PROGRAM  FOR  ACTIVITY  DETERMINATIONS 
i  THE  USAFSAM  WHOLE-BODY  COUNTER 


TED  D  RUPP,  CAPTAIN,  USAF,  BSC 
ROBERT  C.  NELSON,  CAPTAIN,  USAF,  BSC 


NATlOr,--4  TtCHNiCAL 
ink:-'. ‘AA 'ION  SERVICE 


USAF  School  of  Aerospace  Medicine 
Aerospace  Medical  Division  (AFSC) 

Brooks  Air  Force  Base,  Texas 

June  1972 

Approved  for  public  release;  distribution  unlimited. 

'i'h 


lAccasiwi  _ _ 1 

CfSTt 

WHITE  RECTIOX  ^ 

DD3 

BiTF  SKTIOR  Q 

s 

£S. 

□ 

. . . 

1  *  . . 

jer  . 

oisiRiaaiioH/AVCirABiiinf  coob 

OISI. 

A 

AVAIL  3r.i! 

.'6f  smiAi 

Qualified  requesters  may  obtain  copies  of  this  report  from  DDC.  Orders  will 
be  expedited  if  placed  through  the  librarian  or  other  person  designated  to  request 
doeuments  from  DDC. 


When  U.  S.  Government  drawings,  specifications,  or  other  data  are  used  for 
any  purpose  other  than  a  definitely  related  Government  procurement  operation,  the 
Government  thereby  incurs  no  responsibility  nor  any  obligation  whatsoever;  and  the 
fact  that  the  government  may  have  formulated,  furnished,  or  in  any  way  supplied 
the  said  drawings,  specifications,  or  ether  data  is  n'>t  to  be  regarded  by  implication 
or  otherwise,  as  in  any  manner  licensing  the  holder  or  any  other  person  or  corpora¬ 
tion,  or  conveying  any  rights  or  permission  to  manufacture,  use,  or  sell  any  patented 
invention  that  may  in  any  ray  be  related  thereto. 


^curity  Classification 

DOCUMENT  CONTROL  DATA  -  R  &  D 

(Security  f^bittact  end  indexing  unnotHtlon  must  be  enforce/  when  the  ovoratt  report  Is  closstiied) 

I  ORIGINA  TING  ACTIVITY  (Corporate  author)  ze,  REPORT  security  classification 

USAF  Sdiool  of  Aerospace  Medicine  Unclassified 

Aerospace  ^'edical  Division  (AFSC)  it.  croo» - 

Brooks  Air  Force  Base,  Texas  78235 

3.  REPORT  TITLE 

CO'IPUTER  PROURAM  FOR  ACrTVITY  DETERMINATIONS  IN  THE  USAFSAM  WHOLE -BODY  COUNTER 


4  DESCRIPTIVE  NOTES  (Typo  of  report  and  Inctusive  datee) 

Final  report  January  -  May  1970 

9  AUTHORIS)  (First  namct  middle  i/iltlat$  last  name) 

Tod  i).  Rupp,  Captain,  USAF,  BSC 
Robert  C.  Nelson,  Captain,  USAF,  BSC 


0.  REPORT  DATE 

June  1972 _ 

Contract  or  grant  no 


(74.  total  no.  of  PAGES  l/b.  NO.  OF  REFS 


94.  ORIGINATOR'S  REPORT  NUMBER(S) 


6.  PROJECT  NO. 


SAf‘*-TR-72-10 


c.  Task  No.  775701 


OTHER  REPORT  NOIS)  (Any  other  numbers  that  may  be  asslunod 
thia  report) 


10  OISTRIOUTION  STATEMENT 

Approved  for  public  release;  distribution 

unlimited. 

n  SUPPLEMENTARY  NOTES 

13  ABSTRACT 

12.  SPONSORING  MILITARY  ACTIVITY 

USAF  School  of  Aerospace  Medicine 

Aerospace  ^fedical  Division  (AFSC) 

Brooks  Air  Force  Base,  Texas  78235 

- i -  ■  .  .  — 

Tljc  USAFSAM  whole-body  counter  requires  very  precise  calibration  if 
it  is  to  determine  with  accuracy  the  amounts  of  specific  radioisotopes 
contained  in  a  human  body.  Rigorous  detection  efficiencies  for  the  isotopes 
of  interest  must  be  established  in  all  channels  of  the  instrument  and 
frequently  monitored.  Calculations  for  accurate  calibration  and  sample 
determination  involve  lengthy  matrix  operations  wliich  are  too  time-consuming 
for  repetitive  manual  procedures.  Presented  is  a  computer  program  for  the 
accurate  computation  of  the  amount  of  radioactive  material,  as  measured  by 
the  USi.FSAM  whole-body  counter,  in  a  human  body.  Ihe  approach  and  vievelop- 
ment  of  tlie  program  and  matliematical  techniques  are  discussed,  and  all 
factors  involved  in  the  program  are  described  in  detail. 


DD  ;r  J473 


Unclassified _ 

brni'ilN  C'l.i<isificitlion 


Security  Classification 


IVilole-body  counter 
Computer  program 
Calibration 

Activity  determinatioi^j 


ROLE  I  WT 


Unclassified _ 

Security  Classification 


COMPUTER  PROGRAM  FOR  ACTIVITY  DETERMINATIONS 
IN  THE  USAFSAM  WHOLE-BODY  COUNTER 


TED  D.  RUPP,  CAPTAIN,  USAF,  BSC 
ROBERT  C.  NELSON,  CAPTAIN,  USAF,  BSC 


Approved  for  public  release;  distribution  unlimited. 


FORSWOHD 


This  work  was  done  in  the  Health  Physics  Branch,  Radiobiology 
Division,  under  task  No.  775701.  The  study  was  accomnlished  between 
January  and  May  1970,  The  paper  was  submitted  for  nublication  on 
iO  March  1972. 


This  report  has  been  reviewed  and  is  approved. 


EVAH  F..  GOLm,  Colonel,  USAF,  MC 


Coomender 


ii  <L 


ABSTRACT 


The  USAFSAM  whole-body  counter  requires  very  precise  calibration  if 
it  is  to  determine  with  accuracy  the  amounts  of  specific  radioisotopes 
contained  in  a  human  body.  Rigorous  detection  efficiencies  for  the  isotopes 
of  interest  must  be  established  in  all  channels  of  the  instrument  and 
frequently  monitored.  Calculations  for  accurate  calibration  and  sample 
determination  involve  lengthy  matrix  operations  which  are  too  time-consuming 
for  repetitive  manual  procedures.  Presented  is  a  conputer  program  for  the 
accurate  computation  of  the  amount  of  radioactive  material,  as  measured  by 
the  USAFSAM  whole-body  counter,  in  a  human  body.  The  approach  and  develop¬ 
ment  of  the  program  and  nathematical  techniques  are  discussed,  and  all 
factors  involved  in  the  program  are  described  in  detail. 
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COMPUTER  PROCR>\M  FOR  ACTIVITY  DETERMINATIONS 
IN  THE  USAFSAM  WHOLE-BODY  COUNTER 


I .  INTRODUCTION 

Constant  vigilance  is  required  to  maintain  an  accurately  calibrated 
whole-body  counter.  Mass  suppression  and  efficiencies  for  the  instrument 
must  be  continually  remeasured  and  recalculated  to  compensate  for 
component  aging  and  electronic  instabilities.  The  calibration  procedure 
used  is  detailed  by  Dayton  (1)  and,  in  a  somewhat  briefer  form,  by 
Taboada  (2).  A  generalization  of  the  procedure  will  be  presented  to 
illustrate  the  data  required  by  the  computer  program.  Tlie  objective  of 
the  program  is  to  simplify  data  manipulation  after  collection,  thus 
making  possible  a  daily  or  weekly  update  of  calibration  data  points. 

II.  CAIN  AND  WINDOW  SETTINGS 

Calibrated  radioactive  sources  were  inserted  in  the  counter.  Gains 
and  windows  for  each  channel  were  adjusted  to  obtain  maximum  efficiency 
for  the  primary  radionuclide  and  minimum  efficiency  for  other  radionuclides 
which  right  be  present.  Many  criteria  may  be  used  to  define  which  window 
settings  are  best,  but  final  determination  of  the  window  settings  can  only 
be  made  by  trial  and  error,  (Further  information  on  gain  and  window 
settings  is  given  in  references  1  and  2.) 

III.  MASS  SUPPRESSION 

The  addition  of  mass  in  the  whole-body  counter  causes  an  alteration 
in  background.  This  alteration  is  a  function  of  background  energy,  the 
Compton  scattering  cross  section,  and  the  mass  of  tmterial  in  the  counter. 

In  the  energy  region  above  0.5  Mev  the  background  count  rate  tends  to 
decrease  as  the  amount  of  material  is  increased,  vdiile  in  the  region  from 
0,1  to  0.5  Mev  little  or  no  change  occurs.  Tliis  effect  is  caused  by  high- 
energy  gamma  rays  interacting  with  the  material  by  Compton  scattering, 
producing  lower  energy  secondary  radiation.  Ihc  magnitude  of  the  effect 
is  a  few  percent  of  the  net  count  rate  from  a  normal  man;  therefore,  it  should 
be  considered  if  accuracy  to  within  a  few  percent  is  required.  The  effect 
has  been  inr;luded  in  the  computer  program. 

Sugar  phantoms  free  of  gamma- emit ting  radionuclides  were  used  to 
obtain  mass  suppression  data,  since  sugar  is  easily  obtained  find  handled 
and  effectively  approximates  the  elemental  composition  of  a  human  being. 

A  background  count  was  taken  before  and  after  each  phantom  count  to 
compensate  for  any  short-range  electronic  instabilities.  Tne  number  of 
times  a  phantom  is  counted  for  any  given  weight  is  left  up  to  the 
discretion  of  the  user.  The  main  interest  is  to  determine  the  regression 
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line  between  mass  suppression  and  phantom  weight  for  each  channel.  The 
method  of  least  squares  (3)  is  used  to  calculate  the  constants  a  and  b 
in  the  equation  y  =  bx  +  a. 

Further  explanation  can  be  more  easily  accomplished  mathematically. 

To  do  this  some  symbols  need  to  be  defined. 

k  =  subscript  denoting  channel  number  being  calibrated 

j  =  subscript  denoting  phantom  weight 

i  =  subscript  denoting  an  individual  phantom  count 

Wji^  =  phantom  weight,  identified  by  clianncl 

^ijk  “  gross  count  rate  of  phantom  count,  identified  by  channel 
^lijk  background  count  rate  before  phantom  count 
^2ijk  “  background  count  rate  after  phantom  count 

nji^  =  number  of  counts  taken  of  each  phantom  weight  in  a  given  channel 

nj^  =  number  of  phantom  weights  for  each  channel 

yjj,  =  average  net  count  rate  for  phantom  weight  (wjjj) 

Tlic  average  not  count  rate  for  Wj|;  is  given  by  equation  1. 


Tlie  slope  bj^  (c.p,m./lb.)  and  intercept  a)^  (c.p.m.)  for  each  channel  k 
are  found  by  fitting  the  points  (wj^.,  yji^).  The  formulas  for  aj^  and  bj. 
are  given  in  equations  2  and  3. 


*) 


(2) 


and 


Ol  = 


(I«fk)(Iyik)  -  (2*jk) 

i=l  i=i  i=i  j=i  _ 


where 


Dk  '  "k 


"k  „  SU  i 

(K)-  &) 

J=1  j=l 


(3) 


(4) 


The  standard  deviations  for  bj^  and  aj^  are  giveu  by  equations  5  and  6. 

■^w,!  .1/2 


^'’k  °  (^yjkHk 


)LIiL^ 


(5) 


and 


1/2 


(6) 


where 


S. 


,2(>'|k'''’k'«|k+“k’)  ,1/2 
- ) 


(7) 
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which  is  the  estimate  of  the  standard  error.  The  slopes  and  intercepts 
along  with  their  standard  deviations  were  stored  in  the  computer  memory 
for  future  use. 


IV.  I-FFICIENCTES 

The  Final  step  in  the  calibration  is  the  data  collection  and  compu¬ 
tation  of  efficiencies  for  each  channel.  There  will  be  n  efficiencies 
for  each  channel-one  efficienc)'  for  the  radionuclide  of  primaiy  interest, 
and  an  efficiency  for  each  of  the  other  radionuclides  for  ^^hich  the  counter 
is  being  calibrated.  In  all,  there  will  be  n^  efficiencies  when  n  channels 
are  being  calibrated. 

The  efficiencies  and  cross  efficiencies  are  a  function  of  the  phantom 
weight,  as  was  the  background  suppression.  The  gamma  rays  emitted  from 
the  calibration  sources  are  absorbed  and/or  degraded  by  Compton  scattering, 
the  photoelectric  effect  and  pair  production  resulting  in  their  complete 
absorption  or  scattering  into  a  lower  energy  window. 

The  data  collection  and  data  reduction  for  efficiency  determination 
were  very  similar  to  those  used  for  background  suppression.  Background 
counts  were  taken  before  and  after  each  set  of  phantoii^.  The  results  were 
averaged  and  a  best-fit  linear  regression  line  was  computed  by  the  method 
of  least  squares.  The  number  of  phantom  counts  was  determined  by  the 
activity  of  the  standards  and  s'atistics  desired.  Tlie  mathematical 
expressions  wore  similar  in  form  to  the  mass  suppression  equations  with 
the  addition  of  one  more  subscript.  The  symbols  used  are  defined  as  follows: 

k  =  subscript  denoting  channel  number  being  calibrated 

I  =  subscript  denoting  radionuclide 

j  =  subscript  denoting  phantom  weight 

i  =  subscript  denoting  an  individual  phantom  count 

Wjik  =  phantom  weight,  identified  by  channel 

X|  =  decay  constant  for  radionuclide 

S^|]  =  activity  of  standards  corrected  for  decay 

Sjii^  =  activity  of  standard  at  calibration  date 

=  efficiency 

tjii^  =  time  from  date  of  calibration  to  date  of  count 

Yijtk  =  gross  coujit  rate  of  phantom 

^liilk  "  background  count  rate  before  phantom  count 
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-  background  count  rate  after  phanton  count 
Njij^  =  nunfcer  of  counts  taken  for  each  phantom  weight  in  a  given  channel 
Nj)^  =  number  of  phantom  weights  for  each  channel 
^jlk  “  average  net  count  rate  for  phantom  weigjit  Wj|jj 
The  average  net  count  rate  for  weight  Wjjj^  is  given  by  equation  8. 


where  bj^Wjij^  +  a^_  is  the  mass  suppression  correction  term  in  the  background. 
The  equation  for  efficiency  is 


Ejik  "  ^Jlk'^^jlk  where  Sj,^  =  3j||^exp{- A|tj|k)  ( 

The  slope  B|j(  (eff./lb.)  and  intercept  A|j.  (eff.)  are  given  by  equations 
10  and  11,  respectively. 


(11) 


where 


(12) 


The  standard  deviations  for  Bjj^  and  A|j^  are  given  by  equations  13  and 
14. 


^^JlklV/jlk 


(13) 


■  ^EjmlWjik 


(14) 


where 


N|k 

V/  n2\I/2 

ZjVEjlk  "  A|k  Wj||(  “  B||^  )  \ 

/ 


(15) 


is  the  estimate  of  the  standard  error. 
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V.  COMPUTATION  OF  SUBJECT  ACTIVITY 


The  basic  equation  for  calculating  the  activity  for  two  channels  and 
two  radionuclides  is  given  by  both  Da3'’ton  and  Taboada  (1.  2).  A  straight¬ 
forward  generalization  gives  equation  16. 


Sj  +  ^Sl  -  (bii  +  bi2)/2 

Ell  Ei2  ...  liik 

$2  +  MS2  -  (b2i  +  b22)/2 

- 

E2I  ^22  •••  ^2k 

•  • 

■^2 

Sk  +  MSk  -  (bki  +  bk2)/^ 

•  • 

_E  1  j  .  •  •  •  •=  1 1  _ 

1 

■  < 

_ J 

where  the  symbols  are  defined  as  follows: 


(16) 


MS]^  =  bk  W  +  ak 

E|k  =  B|k  W  ^  A,k 

IV  =  subject  weight 

Sk  =  subject  gross  count  rate 

bkj^  =  background  count  rate  before  subject  count 

hk2  =  background  count  rate  after  subject  count 

Ak  =  activity  present 


Tlie  unknown  in  the  equation  is  Ak.  The  solution  can  be  found  by  using 
Cramer's  rule  (4).  The  resulting  determinants  are  evaluated  by  the  triangle 
method  (4). 


Even  though  the  solution  for  Ak  is  straightforward,  it  requires  the 
evaluation  of  a  large  number  of  deteiminants — an  operation  which  is 
tedious  and  prone  to  error.  Tlie  calculations  are  further  complicated  by 
the  error  analysis.  Each  time  an  operation  in  the  evaluation  of  Ak  is 
performed,  the  corresponding  standard  deviation  is  computed.  TTiis  is 
done  until  the  value  and  standard  deviation  for  Aj.  are  obtained.  Three 
simple  equations  were  used  in  approximating  the  standard  deviations: 


2  1/2 

S  =  (si“  +  S2  )  — for  addition  and  subtraction 

S  =  (xi^S2^  +  — for  multiplication 


and 


c  - 


(Cr!)' 


2  1/: 


^  _^for  tho  division  (xj/xt) 


(17) 

(18) 


(19) 
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VI.  COMPUTER  PROGRAM 

The  computer  pro^^-am  consists  of  the  mainline  program  which  calculates 
mass  suppression,  efficiencies,  and  the  table  of  normalized  efficiencies 
versus  weight,  Tlie  mainline  program  also  controls  access  to  the  subroutines. 
Four  subroutines  were  established;  (1)  the  subroutine  BSSS,  which  calculates 
the  activity  and  associated  errors  of  the  subject;  (2)  the  subroutine 
iJECAY,  which  decays  the  standards  from  time  of  calibration  to  present 
date;  (3)  the  subroutine  WLSCF,  which  performs  tlie  least-square  curve 
fit  on  the  mass  suppression  and  effi'^iency  data;  and  (4)  the  subroutine 
EUET,  which  takes  determinants  required  in  the  other  parts  of  tlie  program. 

A  list  of  FORTRAN  symbols,  their  mathematical  equivalent,  if  any,  and 
an  explanation  of  their  use  are  given  in  tables  I  througii  V.  The  same 
symbol  may  be  used  in  several  parts  of  the  program.  A  complete  FORTRAN 
listing  of  the  program  is  in  the  appendix. 

Tlie  data  cards  for  the  program  are  composed  of  eleven  types.  The 
information  contained  on  each  card  t)’pe  and  the  format  for  the  types  are 
given  in  table  VI.  The  program  deck  is  organized  as  shovs-n  in  figure  1. 

Figure  2  illustrates  the  arrangement  of  a  typical  data  deck  for  two  cliannels. 

Tlie  output  of  the  program  will  vary  depending  upon  the  value  used 

for  ITAB.  The  two  main  items  of  interest,  however,  are  the  table  of 

efficiencies  and  the  subject  activity. 

TTie  format  of  the  efficiency  table  is  controlled  by  card  types  9  and 
10,  The  variables  are  always  piinted  in  the  same  order.  The  weight  is 
given  first,  then  the  mass  suppression,  the  efficiencies,  and  the  deter¬ 
minant  of  the  efficiency  matrix.  The  order  of  the  printout  corresponds 
to  the  sequence  of  input  of  the  data  to  the  computer. 

The  subject  printout  is  fixed.  An  example  is  shown  in  figure  3. 

Tlie  program  was  used  to  compute  the  activity  of  approximately  one 
thousand  subjects.  The  results  of  the  program  are  very  reliable  and 
consistent  with  other  methods  of  lean  body  mass  and  potassium  determination. 
TTie  errors  for  the  subject  activity  are  sensitive  to  any  mistakes  in  the 

input  data  and  appear  to  be  a  good  measure  of  the  accuracy  of  the  subject 

activity. 
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TABLE  I 

Symbol  table  for  mainline  progran 


FORTRAN 

symbol 

Mathematical 

equivalent 

Definition 

ABACK (LS) 

"ajk  >>111  Ik 

Background  count  rate  before 
phantom,  divided  by  100 

AIN(I,L) 

^Ik 

Efficiency  regression  line  intercept 

AINT(I) 

Mass  suppression  regression  line 
intercept 

BACK(LS) 

’’2ijk  ^2ijlk 

Background  count  rate  after 
phantom,  divided  by  100 

D 

yjk  0^  yjik 

Net  count  rate  for  phantom 

DAY (IK) 

Julian  date  of  standard  calibration 

DSS(M) 

Determinants  of  efficiencies 

EDSS(M) 

Standard  deviation  in  DSS(N0 

ERRI(I) 

s_, 

ak 

Standard  deviation  in  AINT(I) 

ERROR I ( I, L) 

^Alk 

Standard  deviation  in  AIN(I,L) 

ERRORS(I,L) 

^Bik 

Standard  deviation  in  SLOPE (I, L) 

ERRS(I) 

Sbk 

Standard  deviation  in  SLOP (I) 

F 

Equal  to  NUMB(J)  or  NUMP(N) 

FORM 

Format  for  efficiency  table 

FORS 

. 

Format  for  errors  in  efficiency 
table 

HD(KN) 

Heading  for  efficiency  table 

ID 

Equal  to  NUMB(J) 

IDUM 

Counter 

IDUV 

Counter 

I  DUS 

Counter 
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FORTRAN 

Symbol 

IDUT 

IF 

IG 

IK 

IR 

IS 

ITAB 

lU 

IX 

J 

K 

KJ 

KSSD 

KZZ 

L 

LAMBDA 

LS 

N 

NN 

NUMB(J) 

NUMBC 


TABLE  I  (contd. ) 
Mathematical 

equivalent  Definition 

Counter 

Counter 

Counter 

Counter 


"jk 


Equal  to  NUMP(N) 

Equal  to  NUMBP(N) 

Print  and  computation  control 
Counter 


Equal  to  NUMPO(I) 

Counter 

Counter 

Counter 

Counter 

Counter 

Counter 

Decay  constant  for  standards 

Coup.ter 

Counter 

Counter 

Number  of  counts  taken  for  a  given 
phantom  weight  in  channel  .k 

Number  of  channels  being  calibrated 
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TABLE  I  (contd.) 


FORTRAN 

Mathematical 

symbol 

equivalent 

Definition 

NUMBP(N) 

"k 

Number  of  phantom  weights  per 
channel 

NUMP(N) 

Njlk 

Number  of  data  points  for  phantom 
weight  WGTS(J) 

NUMPO(IK) 

Number  of  phantom  weights  used 
to  determine  efficiency 

PDAY(IK) 

Present  Julian  date 

POINCLS) 

Gross  count  rate  of  phantom 
divided  by  100 

PYEAR(IK) 

Last  two  digits  of  present  year 

SKRCIF,IG) 

Efficiency  matrix  for  weight  SW 

SKT(IF,IG) 

Error  matrix  for  weight  SW 

SLOP(N) 

bk 

Slope  of  mass  suppression 
regression  lines 

SLOPE (N,L) 

®lk 

Slope  of  efficiency  regression  line 

STD (IK) 

®jlk 

Standard  at  time  of  calibration 

STDD(J) 

^jlk 

Standard  corrected  for  decay 

SW 

Weight  in  pounds  for  TABLE (I, J) 

TABLE ( I, J) 

Table  of  efficiencies  and  mass 
suppression  as  a  function  of  weight 

WGT(J) 

'^jk 

Phantom  weights  used  in  mass 
suppression 

WGTS(J) 

W 

.  jlk 

Phantom  weights  used  in  efficiencies 

YEAR (IK) 

Last  two  digits  of  year  source  was 
calibrated 

YS(J) 

Average  net  count  rate  for  phantom 
weight  WGT(J) 

YSS(J) 

yjik 

Average  net  count  rate  for  phantom 
weight  WG’"S(J: 
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TABLE  II 

Symbol  table  for  subroutine  BSSS 


FORTRAN  ^fathe^natical 

synibol  equivalent  .  Definition 


ABACK  (I)  bj^2  Background  count  rate  after 

subject,  divided  bv  100 


BACK  (I)  bj^j  Background  count  rate  before 

subject,  divided  by  100 


CW(K) 

DSS(M3 

DW(K) 

EB(IL) 

EDSS(N0 

ES(IL) 

IJ 

IL 

IM 

ISUB 

IWCT 

J 

K 

M 

N 

NCHAN(I)  k 

NODE 

NUDE 


Standard  deviation  error  in  XS(K) 

Determinant  of  efficiencies 

Standard  deviation  in  SCK) 

Standard  deviation  in  average 
background 

Standard  deviation  in  DSS(M) 

Standard  deviation  in  SOURCE(I) 

Counter 

Counter 

Counter 


IWGr-70 

Subject's  weight  rounded  to  the 
nearest  pound. 

Counter 

Counter 

Counter 

Counter 

Channel  number 

Designates  if  subject's  weight  is  in 
pounds  or  kilograms 

Last  cr'd  indicator 
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TABLE  II  (contd.) 


FORTRAN  MfttheEBRtical 

smbol  eflui  valent 

'i^t^3f*tnntMtr0m>m . .  '■  »»"»«w>iir«afMMWiiMmwiimni— fKniiwM' 


NIIM8C  n 

f.njMS 

S(K)  A|£ 

SOURC(I)  Sj^ 

SOURCE (I) 

TABLE (ISUB»  IJ) 


TIME 

WCT  W 


XCN,M)  Eju 

XS(K) 

XX(N,M) 

Y(N,M) 

Z(K,N) 


Definition 


Number  of  channels  being  calculated 

Subject  sample  number 

Disintegrations  per  minute, 
divided  by  100,  for  channel  k 

Gross  count  rate  of  subject, 
divided  by  100 

Net  count  rate  of  subject 

Table  of  mass  suppression 
efficiencies  and  their  standard 
deviations 

Length  of  time  for  count 

Subject's  weight  in  pounds  or 
kilograms 

Efficiency  matrix  for  subject 

Determinant  of  Z(K,N) 

Working  storage 

Standard  deviation  matrix  for 
X(N,M) 

Working  storag*. 


& 


; 

I 

r 


TABLE  III 

Symbol  table  for  subroutine  DFiCAY 


FORTRAN 

synbol 

Mathematical 

equivalent 

Definition 

CCURIE 

Sjik 

Source  strength  corrected  for 
decay 

CURIE 

Sjik 

Source  strength  at  time  of 
calibration 

DAY 

Julian  date  source  was  calibrated 

LAMBDA 

^1 

Decay  constant 

PDAY 

Present  Julian  date 

PYEAR 

Last  two  digits  of  present  year 

XPDAY 

Number  of  days  between  PDAY  and  DAY 

XPYEAR 

^jlk 

Elapsed  time  between  calibration 
date  and  present  date 

YEAR 

Last  two  digits  of  year  source 
was  calibrated 
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V 


TABLE  IV 


Symbol  table 

for  subroutine  WLSCF 

FORTRAN 

Mathematical 

symbol 

equivalent 

Definition 

AINTER 

b 

Y  intercept  of  the  regression  line 

DELTA 

Working  storage 

DIFF 

Difference  between  experimental 
and  calculated  value  of  Y(I) 

ERRORI 

Sb 

Standard  deviation  of  AINTER 

ERRORS 

Sa 

Standard  deviations  in  SLOPE 

I 

Counter 

N 

Counter 

SLOPE 

a 

Slope  of  the  regression  line 

SN 

Working  storage 

SUMX 

lOy. 

E  Xi 

i«l 

10 

Sum  of  X(I)’s 

SUMXY 

Z  XiYi 
i-l 

10 

Sum  of  the  product  X(I)  and  Y(I) 

SUMY 

?  Ti 

Sum  of  Y(I)  's 

i»l 

SY 

X(I) 

•  Xi 

Independent  variable 

Y(I) 

Yi 

Dependent  variable 

TABLE  V 

Sypfcol  table  for  subroutine  EDET 


FORTRAN 

syrfjol 

Mathemtical 

equivalent 

Definition 

B(KS,K) 

Reduced  matrix 

C 

aij/ajj 

Reducing  multiplier 

D 

Det|aij| 

Determinant  of  S(KS,K) 

ERB(KS,K) 

Standard  deviation  in  B(KS,K) 

ERB2 

Variance  in  B(KS,K) 

ERD 

Sd 

Standard  deviation  in  D 

ERD2 

Working  storage 

ERC2 

Standard  deviation  in  C 

ES(KS,K) 

®aij 

Standard  deviation  matrix  for 
SCKS,K) 

J 

Order  of  matrix  S(KS,K) 

JK 

Counter 

JS 

Counter 

K 

» 

Counter 

KI 

Counter 

KR 

Counter 

KS 

Counter 

KT 

■ 

Counter 

L 

Counter 

S(KS,K) 

Matrix  to  be  evaluated 
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Innut  data  cards 


Type 

Columns 

1 

1-50 

51-80 

2 

1-3 

4-6 

7-9 

10-12 

13-15 

16-18 

19-21  ' 
22-80 

3 

1-4  ' 

5-7 

8-11 

12-14 

J.5-18 

Variables 


1-9 

10-16 

17-23 

24-30 

31-70 

71-80 


Foritat 


50H 


NITMBC 

NU^fP(l) 

NU^^P(2) 

NUNfP(3) 

NUMP(4) 

NWfP(5) 

ITAB(*) 


WCT(l) 

NU>fB(l) 

wen  (2) 

NU^fl^(2) 

WGT(3) 


POINCLS) 
BACKCLS) 
ABACK (LS) 


Remarks 


Title 

Unused 


NUMBC^S 


NUMBP(I)<10 


Unused 

Unused 


NUWh'D^IOO 


Repeats  until  subscrip 
reaches  NUMBP(I) 


TABLE  VI  (contd. ) 


Columns 

Variables 

Format 

1-3 

NU^(P0C1) 

13 

4-6 

NUMl’0(2) 

13 

7-9 

NUMP0(3) 

13 

10-12 

NlMf>0t4) 

13 

\ 

13-15 

NUMPOCS) 

13 

16-22 

LAMBDA 

F7.6 

23-80 

1-4 

WGTS(l) 

F4.0 

5-11 

STD(l) 

F7,2 

12-13 

YEAR(l) 

F2.0 

14-17 

DAY(l) 

F4.0 

18-20 

NUMP(l) 

13 

21-24 

WGTS(2) 

F4.0 

25-31 

STD (2) 

F7.2 

32-33 

YEAR(2) 

F2.0 

34-37 

DAY (2) 

F4.0 

38-40 

NUMP(3) 

13 

1-9 

9X 

10-16 

POIN(LS) 

F7.1 

17-23 

BACK(LS) 

F7.1 

24-30 

ABACK (LS) 

F7.1 

31-48 

18X 

49-51 

PYEAR(LS) 

F3.0 

52-54 

3X 

- -  -  « 

Remarks 


NUMP0'(I)<10 


Floating  point 
Unused 


NUMP(I)<100 


Repeats  until  subscript 
reaches  NUMPO(I) 


Unused 


Unused 


TABLE  VI  (contd. ) 


Type 

Columns 

Variables 

Format 

Remarks 

7 

55-57 

PDAYCLS) 

F3.0 

58-80 

> 

Unused 

8 

1-80 

/ 

13A6,A2 

Title  card  for 
efficiency  table 

9 

1-60 

FORM  or  FORS*^ 

10A6 

Format  for  printout 
of  efficiency  table 

10 

1-78 

f 

78H 

Subject  name, 

SSAN,  date,  etc. 

11 

1-4 

NU>fS  / 

14 

S-12 

WGT 

F8.5 

13 

NODE 

11 

14-22 

TIME 

F9  4 

23 

NCIIAN(l) 

11 

NC»AN(I)^NUMBC 

24-30 

SOURCE  Cl) 

F7.1 

31-37 

BACK'Cl) 

F7.1 

38-44 

ABACK(l) 

F7.1 

45 

NCUAN(2) 

11 

Repeats  until  subscript 
reaches  NIWBC 

*ITAB  controls  the  function  of  the  program  according  to  the  following 
conditions; 


M3DCITAB,2)5  0 
mU(ITAB,4)  ?2 
M0D(ITAB,16)^8 

M0D(ITAB,32)?16 
WDCITAB,64)«2' 
/input  cards  not  required 
/input  cards  not  required 


Print  table  of  efficiencies  as  a  function  of  weight 
Calculate  activity  of  subjects 
Print  slopes  and  intercepts  for  mass  suppression 
and  standards 

Print  mass  suppression  data 
Print  standard  calibration  data 
unless  MOD(ITAB,2)50 
unless  MOD(ITAB,4)3;2 
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'YPE  3-COWTROL 
PE  2 -MASS  SUPI 
E  1-TfTLE  CARD 


DOE,  JOHN  E.  OOOOOOOOO(^) 

1  JAN  70 

0000 200.0^®^  CHANNEL  NO. 

DPH/100 

ERROR/100 

1 

186.85^^^ 

18.92^^^ 

2 

268.02 

27.24 

(®)  Subject’s  name. 

Social  security  number. 

Date  count  was  made. 

Subject’s  code  number. 

(®)  Weight  in  pounds. 

Activity  divided  by  100.  Tlie  units  on  the  activity  are  the 
same  as  those  on  tlie  standards.  They  must  be  consistent  throughout 
the  program. 

Standard  deviation  in  the  activity,  divided  by  100, 

FIGURE  3 

Example  of  subject  output. 
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APPENDIX 

FORTRAN  LISTING  OF 
WHOLE-BODY  CDUNTER  PROGRAM 


Preceding  page  blank 


25 


*  ^  in 

N 

in 

•  o  a 

in  I  O 
—  •  a 
z  - 1: 

1  m  uj 

I  •  tn 

;  ^  —  in 
:  m  f-  •> 

I  •  m 

■  m  i/)  — 

.  —  •  i/) 

I  llI  a: 

■  a  m  o 

.  O  •  Qi 
I  j  in  O' 

I  l/)  w  ui 
)  •■  O'  - 
>  ~ 

•  O  tO  •-< 

>  t-«  *  o 
4  w  ^ 

.  J)  o  fi 

rg 

C  o 

j  3:  -  - 

•  *  >-  UJ 

►  -j 

'  o  o  00 

)  rH  Q.  < 
H  •  H- 

*  Q.  —  • 

0  Z  O  — 
7  3  0'-* 

-  z  r-(  nj 
»  »  nj 
^  —  q;  — 
3  in  <  i/) 

H  —  UJ  U) 

-  O  >-  O 

-  0.  a  •> 

I  s:  .  ~ 
3  3  -*  <H 
•>  Z  O  nj 
►  iH  nj 
3 - - 

-4  m  lO  u) 

-  —  q:  ui 
/)  Q  O  O 
M  t-  U.  UJ 

*  to  ► 

n  o  o 

in  *H  tH 

X 

D  2  O 
E  q:  a  Q 

D  OJ  o  »- 

Z  UJ  U.  (/) 
-4  (M  (<><}■ 


< 

Q  O 

00  o  o 
S  o  o 

<  O  rH 
-J  r-4 

-J  o  z 

<  <  -< 
uj'uj  oi: 
O'  q:  a 


I  o  o  o  o 

I  •  •  •  • 

t  o  o  o  o 

M  II  11  II 
■  Z - - 

>  w  z  z  z 

I  Q.  w  w  w 

« CO  a  to  »- 
z  o  q:  z 
)djq:-* 

)  Z  lOUJ  < 


O  O  O  r-4 

•  •  •  II 

000 

II  11  II  z 

~  z 
z 

z  ^  w,  <n 

—  zoo 

—  V  Q.  iH 

0:  o  z: 
q:  I-  o  o 
LU  lO  z  o 


o 

o  •  o 

•  o  • 

O  II  o 
II  II 


z  z  z 
z  -  z 
q:  — 
to  <  >- 
lO  UJ  < 
Q  >-  Q 

UI  a  a 


o 

o  o  O  f-i 

•  •  •  » 

O  O  O  f-i  o 

II  II  II  II  • 

—  o 

- Z  Z  II 

z  z  z 
z  z  —  m  — 
- o  z 

Z  5i^  VJ  »H  — 

•-I  u  <  «o 
o  <  CD  O  10 
a  CO  <  o  >■ 


u  u  u 


Preceding  page  blank 


f 


o  o 

•O  *00000 

o»o»»»#* 

II  O  II  o  o  o  o  o 
II  II  II  II  II  II 

z  z  "«  z  z  z 

-Z'-Z'- - - 

Q'-a  —  0O2Q.I/) 

>-<UJOI5</)3e) 

iOQ>-3Z>-Z3: 

m 

o 


Q 

Q  -. 


(/) 

z 

o 

1-^ 

t- 

< 

_l 

D 

U 

_l 

< 

u 

z  • 

o 

•—I  • 

10  • 
UJ 

cn.  • 
a 

a  • 

D 

lO  • 

If)  • 

»0 

<  • 

z 


CO 

•  < 


•  ~  u 

U  CQ 

•  CO  z 

Z  D 

•  Z 

z  o 

•  »  \0  f-t 

II 

•  II  o  ►-* 

>-*  h—  •» 

•  * 

^  o 

■  r-l  W 

Q_ 

•  a  —  M 

CO  o  z 

•  Z  t-i  3 

o  •  z 

•  z  I-  — 


II 

10 


II 

10 


u  u 

CO  00 

z  z 
z  z 


u  — 

CO  (\J 

z  m 
o  •> 
z  CO 

iH  »-« 
rH  f—  O  rH  w 


o  -< 

O 

rH  Q 

O 
Q  Z 
< 

UJ  u. 
oc  *-. 


O  II 


t—  rsj 

Z  -t 


o 

Q 

O 

vO 


lO  •-• 

►-i  rH  •» 

•  \0  rH 

•H  II 

HOT 

-1  h-  •> 

~  O  “> 

“)  (5  w 

—  m 

CO  ~  z 

Z  O  D 

D  rH  Z 

z  •  • 

•  I - - 

~  -J  “) 

->  •  w 

M  r-^  ^ 

I-  t\J  o 

o  cn  X 

3  •  w  lO 

—  OO  •  >-H 

•  <  1\J  •> 
(M  I-  O  rH 

O  -I  o  II 
O  rH  -) 

rH  Q 

O  I-  O 
Q  Z  Z  rH 

<  —  H-r 

UJ  u.  or  o 
q;  r-,  0.  o 


10 

_j 

v; 

u 

< 


lO 
_j 

u 
< 
CO 

CO  CM  < 
<  \0  » 
•  r» 

~  O  lO 

U)  I-  _J 

_l  - 

—  o 
o  u 
u  < 

<  —  CO 
CO  xo  • 

»  rH  r. 

r*  •  10 
lO  I-  _J 
_J  _J  w 
W  ,  Z 

z  — 


o 

a 


“)  CO 
w  o 
CO  o 
Z  rH 

Z  Q 
II  < 

a  UJ 

*-t  q: 


CNJ  o 
CO  a 

»  w 

CO 

< 


O' 

a 

UJ 


CO  o 

o  • 

o  o 

rH  II 


z 


«  r^  -) 

rH  ^  (/) 


o 

• 

o 

w 

D 

1-  a 

r-H 

z  a 

* 

*-<  UJ 

<  » 

o 

• 

i-H 

CVJ 

w 

>-  CVJ  1- 

lO  r-c  z 

a  •-* 

a  o  < 

*rr 

a  h*  •> 

u 

o  - 

< 

►H  lO  *r 

CQ 

—  10 

< 

a  —  a 

+ 

o  00  a 

a  a 

lO  •  • 

r-r 

h“ 

to  a  rH 

u 

>-  •  >- 

<  Q 

••  r.  Q. 

JO  + 

u. 

1-  vO  o 

N. 

O  rH  _| 

1  ~ 

jE  *>  10 

“)  >- 


vO 


a  a  to 
•-•  a  >- 

rvj 

vO 


II  Z  > 

iiiC  II 

O 
rH  a 

rH  V-  -> 

I  — 
O  II  to 
.  Q  Q  >- 


If) 

>■  U  ■ 

II 

4 

~)  \- 
—  z 
to  o 

>-  VJ 


a 
vj  t- 
tO 

_j  >- 
■S  Q 
O 
-J  Z 


c\j 

cvj 


<  a  a 
u  -<  a 


a  II 
3  — 

z 

o 

CVJ 
Z 

o  o 

U  Q 

CVJ 


Z  H- 


U  U  U 


28 


TABLE(Itl)  =SW 
DO  20  N=1»NUMBC 


• 

• 

* 

•k 

r-4 

• 

II 

II 

to 

to 

• 

_J 

_i 

. 

• 

• 

lO 

to 

• 

X  X 

_l 

-1 

w 

w 

• 

>- 

>- 

< 

< 

• 

II  11 

C3 

a 

r 

r 

a 

a 

• 

. 

• 

~  -> 

lO 

to 

~>  — 

_J 

_i 

fM 

• 

—  a 

* 

a  z 

a 

a 

♦ 

• 

Z  3 

< 

< 

z>  z 

Ul 

UJ 

Z 

• 

z  . 

>- 

>- 

♦ 

a 

a 

• 

"3 

. 

O' 

< 

“)  w 

ar 

• 

a 

>-  >- 

to 

lO 

LU 

00 

>-  < 

_j 

_j 

+ 

• 

5:  < 

<  0 

c\) 

<  Q 

Q  •> 

* 

• 

-J  OQ  O 

U 

u 

-<  ♦ 

•  z  • 

< 

< 

z 

• 

—  <  o 

■3  — 

U3 

J3 

—  36 

00  _l  II 

—  O'  a 

< 

"t  < 

K-  if) 

• 

>0  • 

OC  'O  < 

•> 

\0  . 

Z  ♦ 

< 

<  UJ 

•— <  ^ 

• 

O  in  Q 

UJ  0  >- 

to 

0  to 

<  z 

in  h-  .  00 

>“  1—  . 

_J 

1-  _l 

+  ^ 

• 

.  f-i  2; 

« 

r-i  0  II  < 

—  0  “3 

0  irC 

3  a; 

• 

II  0  _j 

-)  0  >- 

u 

10  u 

U)  QC 

^  t-i 

w  Q 

< 

< 

♦  UJ 

• 

»-4  ^  ^ 

Q  —  1- 

CO 

-» CO 

♦  eg  0 

I-  eg  1/3 

. 

CVJ  . 

z 

• 

—  m  rH 

mm.' 

(<3  — 

W  H- 

.  •-<  0 

to 

.  to 

a  on 

• 

»-<  h-  0 

—  1-  -> 

.  _J 

1—  _J 

o  o 

^  _l  0  0 

~3  -J 

0 

_J  — 

-J  uO 

•  lU 

0  •  a  0 

^  .  1/3 

0  z 

.  z 

1/)  II 

a  '-•  z  u  . 

lO  1- 

1-4  ►— 

^  -H 

II  ~ 

• 

Z  4-  D  33  0  0 

I- 

♦  0 

0 

fH 

33'OZZ  .~r-«o 

0  «o  3: 

—  a 

0  a  a 

+ 

• 

Z*  —  DH-  ►-<.» 

36  .  w  X 

~)  — 

^  *—4 

z  z 

-'00  •Z-J'-'-iO 

W  CQ  »  >-• 

w  . 

CD  .  0  . 

*  ♦ 

• 

.  <  •  .  0  it  II 

•  <  in  • 

0  “3  eo 

<00  .  —  r-4 

rg  (M 

^  h-  0  «-<  <  a 

in  1—  0  1—  »-  0 

1-0  0  ~3  II 

.  . 

• 

0*-*0  II  QZ-*  — 

0  -•  0  II 

lO  a  0 

‘-*0  II 

*-4  t—t 

0— 

0  1-4  ”3 

1!  Z  •-• 

*'*  >-4  a  ^ 

• 

Q  Z  Z 

r-4  0 

—  =3 

Q  z 

UJ  UJ 

0  1—  0  <t  cn  " 

0  H-  in  -3  z 

0  t-  “3  3)  'O 

-J  _J 

• 

ozztn_j  iif-Q 

0  Z  Z  -I 

W  II  0 

Z  Z  ^  z 

(0  03 

0  <  ~ 

W  ^  lO  tl 

<  <£ 

• 

II  ujaaoaxoi- 

ujaaoi-aujaaaaio  0 

f- 

-jQ:'-^ao»-'*-»Qu)a*-«ao 

1/)  a  1-4 

l-l  Q.  4-4  >.  U.  0 

• 

rH 

rH 

O 

rH  GO  CO 

O' 

sT 

r\j 

• 

• 

rg  ^0  r- 

gj 

'O 

u  u  u 


29 


POIN(K)=POIN(K )*100. 

BACK  (  K. )  =BACK  (K  )  *100. 

ABACK(K)=ABACK (K) *100. 

CALL  DECAY { STD ( J ) » YEAR( J ) »DAY( J) tPYEARiK ) *PDAY{K ) » LAMBDA »STDD ( J ) ) 
D=(P0IN{K)-(BACK{K)+ABACK(K) J/2.0+SL0P( I )*WGTS( J)+AINT( i n/STDD( J) 
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